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1. 

CHAPTER  Xm 
INTRODUCTION^ 

1 .Object  of  the  Thesis % 

There  are  a  large  number  of  methods  for  determining  to  any 
desired  degree  of  accuracy  the  complex  roots  of  numerical  equations, 
The  problem  has  been  attacked  by  many  of  the  foremost  mathematicians 
and  many  methods  of  approximation  have  been  developed .This  thesis 
alms  to  present  and  to  discuss  the  more  Important  of  these  methods, 
Btatlng  the  relations  between  the  various  ones  wherever  possible. 

This  work  will  deal  only  with  equations  whose  coefficients 
are  real .Though  some  of  the  methods  treated  can  be  applied  directly 
to  equations  with  complex  coefficients, it  is  not  the  purpose  of 
this  article  to  discuss  this  application  in  detail. The  problem  for 
equations  with  complex  coefficients  is  not  more  general, but  leads 
to  equations  of  higher  degree  with  real  coefficients, it  being 
merely  necessary  to  multiply  the  original  function  by  the  conjugate 
function  In  order  to  make  all  coefficients  real. 

No  attempt  will  be  made  to  go  carefully  into  the  graphical 
methods  for  determining  complex  roots. The  algebraic  methods  will  be 
adhered  to  except  in  those  cases  in  which  a  graphical  representa- 
tion will  make  clear  the  algebraic  processes  in  use. 

An  outline  of  the  study  of  complex  roots  has  been  given  by 
F.CaJori  in  an  article  entitled  WA  History  of  the  Arithmetical 
Methods  of  Approximation  to  the  Roots  of  Numerical  Equations'* 
published  in  the  Colorado  College  Publication  of  November, 1 910. 
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2. 

The  article  deals  primarily  with  the  history  of  the  study  of  all 
roots, but  those  parts  of  it  regarding  complex  roots  have  been 
of  value  in  the  preparation  of  this  thesis. The  references  contained 
in  it  have  been  especially  valuable. 

A  list  of  papers  dealing  with  this  subject  or  with  phases  of 
it  is  given  at  the  end  of  this  thesis. 

2. Historical  Outline,  ir'b'y'f^x 

The  history  of  the  development  of  methods  for  the  determina- 
tion of  complex  roots  of  algebraic  equations  Is  so  closely  bound 
up  with  that  of  real  roots  that  the  separation  of  the  two  is 
rather  difficult .Many  of  the  methods  which  apply  to  complex  roots 
were  discovered  in  the  search  for  a  means  of  determining  real 
roots, though  some  of  them, devised  originally  for  real  roots, were 
not  applied  to  complex  roots  until  many  years  after  their 
discovery, as, for  example, the  methods  of  Newton^    and  of  Horner/3 

The  researches  of  Budan',b   and  Fourier av^in  the  first 
quarter  of  the  nineteenth  century, and  later  the  theorems  of 
Sturm      fo  in  1835, of  Cauchy  "',2>  '3  and  of  Laguerre  39jb  laid  a 
foundation  for  work  on  complex  roots .The  method  of  Newton, first 
developed  in  1671, was  applied  to  complex  roots  by  Simpson^ 67  about 
1750  and  by  Cauchy 0  in  1849.*  The  method  of  Taylor's  theorem, 
proposed7for  real  roots  in  17 17, was  not  applied  to  complex  roots 
before  the  ninteenth  century. 

In  1728  Daniel  Bernoulli #  suggested  the  use  of  recurrent 

series  for  finding  roots  of  equations  and  this  formed  a  basis  for 

*waring  is  credited  by  Cajori,  6  p,  ill,  with  having  first  extended 
Newton's  method  to  complex  roots, but  no  reference  is  given. 
frThls  reference  was  not  avail abi See  Cajori. V- 
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work  don©  later  by  EuleraH  and  Lagrange 3       in  developing  what  is 
now  called  Bernoulli's  method, and, of  still  greater  Importance, for 
the  method  worked  out  by  Graf f e  3i     in  1 837 . 

Since  the  first  part  of  the  ninteenth  century  there  has  been 
a  great  interest  in  the  determination  of  complex  roots  which  has 
resulted  in  the  developement  of  a  large  number  of  new  methods  as 
well  as  the  extension  to  complex  roots  of  methods  for  finding 
real  roots. In  1870  Schroder       made  an  attempt  to  unite  a  number 
of  these  into  one  general  method,His  results, while  perhaps  of 
little  practical  value, are  useful  theoretically  in  that  they  show 
the  relationship  between  several  methods ;while  the  article  is  at 
least  a  step  towards  the  generalization  of  methods  of  attacking 
the  problem, 

3 .General  Classification, 

The  methods  of  approximation  may  be  divided  into  three 
classes, as  follows: 

1 .  Methods  of  repeated  limits(with  rectangular  coordinates). 

2.  Methods  of  repeated  limits(with  polar  coordinates), 

3.  Methods  of  double  limits. 

In  the  first  class  the  root  will  be  expressed  in  the  form 
^~K/3.lt  will  be  found  by  means  of  two  distinct  limiting  processes; 
the    one    will  determine  the  real  part  of  the  root  <t  ,the  other 
the  imaginary  part  <-^3  •  No  matter  whether  ol  or  ^3  is  to  be  first 
determined, it  will  be  definitely  fixed  before  any  step  is  taken 
to  find  the  other, except  in  some  cases  for  a  rough  preliminary 
location  of  the  complex  roots  which  precedes  the  more  accurate 


determination ♦This  method  of  repeated  limits  may  be  expressed  as: 

L    L    (  oi'^Vi')  ■    </-K/3  . 

where  o(-f^3  is  an  approximation  to  the  root 

The  second  class  is  similar  to  the  first, except  that  polar 
coordinates  are  used; thus: 

Alt'  onL6 

or  L     L  f*  «  . 

where   f£    is  an  approximation  to  the  root  • 

The  third  class  contains  those  methods  making  use  of  double 
limits; that  is, neither       nor  /3     (nor  /  nor  0  )  is  definitely 
determined  before  the  other, but  successive  approximations  give 
improved  values  for  both  c(  and  S  (or  ^  and  #  );thus: 

or  L   7/  « 

This  method  of  classification  is  very  convenient  for  the 
purposes  of  this  paper, as  it  includes  all  cases  that  will  be 
consldered.There  are  cases  in  which  the  complex  number  is  not 
broken  into  two  separate  parts, such  as  certain  graphical  methods 
or  methods  of  vectors. These  would  of  course  not  come  under  this 
classificatlon.lt  is  the  object  of  this  thesis, how ever, to  deal 
only  with  algebraic  methods, and  all  of  these  fall  very  conven- 
iently Into  the  classes  as  stated  above. 
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CHAPTER  II. 

METHODS  OF  REPEATED  LIMITS.  (  RECTANGULAR  COORDINATES  ). 

p 

4. Lagrange's  Method.     '  > 
Given  an  equation 

vm,  fd)  «  o, 

of  the  nth  degree, with  real  coefficients, having  two  or  more 
complex  roots      V3j,  U~  VO-  Form  the  equation  whose  roots  are 
dqual  to  the  squares  of  the  differences  of  pairs  of  roots  of 
the  given  equation. The  difference  of  two  conjugate  roots  is 

2. 

21^3  which, when  squared, gives  -4p  .The  square  of  the  difference 
of  two  real  roots  will  be  a  real, positive  number .The  square 
of  the  difference  of  a  real  and  a  complex  number, or  of  two  non- 
conjugate  complex  numbers  will  in  general  give  a  complex  number* 
Hence  to  each  real, negative  root  w  of  the  new  equation  will 
correspond  two  conjugate, complex  roots  of  the  original  equation. 
Also, since  w  ■  -4/3^, the  imaginary  parts  of  these  two  roots  will 

be  equal  to    -lfrCw  ,_i£-w  . 

2  2 

Having  determined  (3    , substitute  c(-h((3   for  z  in  f(z)  m  0, 
obtaining  two  equations  in 

u(*)  =  0,  v(o<  )  =  0, 
where  u  includes  all  the  terms  of  t(<ii~^(3  )  that  are  real, while 
v  includes  only  those  involving  l.The  correct  value  of  that 
corresponds  to  (3  is  obtained  by  finding  the  common  root  of 
u  =  0  and  v  ■  0.  This  can  be  done  by  means  of  rational  operations, 
it  being  equivalent  to  finding  the  greatest  common  divisor  of 


the  functions  u  and  v. 

There  are  two  cases  which  require  special  attention. First, 
if  a  real  root  is  equal  to  the  real  part  of  a  complex  root, the 
difference  between  the  two  will  be  the  imaginary  part  of  the 
complex  root; and  the  square  of  this  difference  will  be  a  nega- 
tive number .Second, if  two  complex  roots  have  equal  real  parts 
but  the  roots  are  non- con Jugate, the  squared  difference  will 
again  be  a  negative  number .In  either  case  we  will  have  a  nega- 
tive root  of  the  equation  of  squared  differences  to  which  two 
conjugate  complex  roots  do  not  correspond. 

To  examine  the  first  case, say  that  the  complex  roots  are 
U-f  i(6)J(<k-i^)J  and  the  real  root  <*  .Then  the  equation  of  squared 
differences  will  have  two     roots  equal  to  -/3X;but  to  Iwu  of 
these  will  not  correspond  two  conjugate  complex  roots  each. 77 

In  the  second  case, let  the  roots  be  {^~l^J-  (d+  lV,),  (o(~ 

Then  the  equation  of  squared  differences  will  have  as  roots 

Here  are  two  pairs  of  equal  roots  to  which  do  not  correspond 
two  conjugate  roots. To  the  single  roots  -4^3  »-4^ , however, there 
are  corresponding  pairs  of  conjugate  roots. 

Lagrange  3  ^7?deduces  the  following  conclusions: ( 1 )When  all 
the  real  negative  roots  of  the  equation  of  the  squares  of  the 
differences  of  the  roots  of  the  original  equation  are  unequal, 
the  original  equation  will  have  a  pair  of  conjugate  complex 
roots  corresponding  to  each  one. (2). If  among  the  negative  roots 
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of  the  squared  differences  equation  equal  roots  are  found, then 

each  unequal  root  will  always  furnish  a  pair  of  imaginary  roots* 

Each  pair  of  equal  roots, however, may  give  two  pairs, or  no, 

imaginary  roots ♦Three  equal  roots, from  similar  reasoning, may 

give  six, or  two; four  may  give  eight , four, or  none; and  so  on*  In 

these  latter  cases, the  exact  number  of  complex  roots  corresponding 

to  equal  negative  roots  of  the  equation  of  squared  differences 

can  be  determined  when  u  =  0  and  v  =  0  are  solved  for   °(  .For 

example, if  we  get  two  equal  roots  w,u  =  0  and  v  ■  0  must  give 

two  equal  values  of  o(  .If  these  two  values  are  real, the  equation 

has  two  pairs  of  complex, conjugate  roots  corresponding  to  the 

two  equal  values  of  w;if  complex, none. 

In  theory, Lagrange  has  developed  a  very  clear  and  exact 

method  for  the  determination  of  the  complex  roots  of  an  equation* 

from  the  practical  side, however, it  is  almost  without  value, at 

least  for  equations  of  higher  degree  than  the  fourth. The  process 

of  building  up  the  equation  of  squared  dif ferences , while  possible 

by  means  of  the  symmetric  functions  of  the  roots, becomes  very 

complicated;and  even  if  one  could  obtain  it  easily, it  would  be 

of  degree  n(n- 1 )    which  for  n  greater  than  4  becomes  so  large 
2 

as  to  render  the  determination  of  it's  real  roots  very  tedious* 

^Example:  Find  the  complex  roots  of 

x**-f  x3-5xL-7x-M0  ■  0. 

The  equation  of  the  squared  differences  is 

x<»-43x^  -679*y  -4209x3  -  248xV  99724x  -  96400^4 

Solving  this  equation  for  real  roots, we  find  them  to  be  -4, and  1* 

We  have  therefore    one  pair  of  complex  roots  of  the  original 
_gfeg_jnr^mt>leg  -given  are  Boivea  pv  several  methods,  where  practical ._ 
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equation  whose  Imaginary  part 

S  =±i£w  gsiitfT  =±1*. 

'  2  2 

Substituting  for  x  in  the  original  equation, and  setting 

p>  =  1  ,we  have 

q  ^  4-2>av-Wo(  -  £   =  0 
The  greatest  common  divisor  of  the  left-hand  sides  is    o(  +  .z  . 
Hence  -2  is  the  root  sought, and  the  two  desired  roots  of  the 
original  equation  are  -2  +  i,  -2  -  i.Had  we  placed  (3  m  -1  in  the 
above  equations, the  resulting  roots  <K  would  have  been  the  same 
which  would  be  due  to  the  fact  that  (?>  enters  to  even  powers 
only  in  the  two  equations . 

5 .Method  based  upon  Taylor's  Theorem,     l°'  i,J  *    7?l  1  ~1' 

Given  an  equation  of  the  nth  degree  with  real  coefficients 

f(z)  m  f(x+  iy)  =  0, 

where  x  and  y  are  real. Expand  by  means  of  Taylor's  theorem , thus : 

f(x4ly)  =  f(x)+  f(x).iy  ,  f"(x)-lzy  +  ...  . 

1!  2!  n! 

Since  the  condition  that  f(x-(~iy)  =  0  is  that  the  sum  of  the 

real  terms  and  the  sum  of  the  imaginary  terms  must  each  be 

equal  to  zero, we  have 

f(x)  -  ffriar3-,  f*»ti>^y  =  o, 

2!      ^  4! 

f '(x)  -  f'MxKv1;  fy(x).y?-  ....  =0, 
1!  3!       +  5! 

where  each  equation  has  a  finite  number  of  terms, and  is  nf 
degree  not  greater  than  n  in  y;  or, since  only  even  powers  of  y 
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appear, of  degree  not  greater  than  n  in  y1, 

8 

Eliminating  y  from  these  two  equations, we  can  solve  the 

resulting  equation  for  real  values  of  x, which  we  can  substitute 

in  one  of  the  above  equations  to  obtain  y.The  equation  in  x  will 

have  n(n-1 ]  roots, each  equal  to  one-half  the  sum  of  a  pair  of 
2 

roots  of  f(z)  =  Of Had  we  eliminated  x,the  resulting  equation  in 
y-2-  would  have  had  roots  equal  to  -J.  the  square  of  the  differences 

of  pairs  of  roots  of  f(z)  =  O.To  determine  roots  of  f(z)  ■  0, 

we  need  consider  only  real  values  of  x  and  y,all  such  pairs  for 

which  y  5*  0  giving  complex  roots. 

This  method  is  evidently  closely  associated  with  that  of 

Lagrange, and  is  subject  to  much  the  Same  limitation. The  process 

of  elimination  between  the  the  two  equations  is  very  complicated, 

leading  as  in  the  former  case  to  a  solution  of  an  equation  of 

degree  n(n- 1 )  ;  in  fact, if  x  is  the  variable  that  is  eliminated, 
2 

the  equations  to  be  solved  have  identical  roots  except  for  a 

constant  factor.  Hence  this  method  has  little, if  any, advantage 

over  that  of  Lagrange  if  judged  from  the  practical  standpoint . 

Where  one  method  demands  the  formation  of  the  equation  of  squared 

differences  and  the  determination  of  the  greatest  common  divisor 

of  two  functions, the  other  requires  the  elimination  of  a  variable 

between  two  equations  and  later  the  solution  of  an  additional 

equation  for  it's  real  roots. 

The  two  equations  in  x  and  y  could  have  been  obtained  by  a 

substitution  of  x+  iy  for  z,and  the  separation  made  as  before»The 
*  See  Bauer, '    p. 239* 
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problem  becomes  the  same  throughout •  ' 

Example:  z^-f-z1  -5zx-7z-f10  ■  0. 

By  Taylor's  Theorem, 

f(x-My)  ■  f(xK  f'(x)»ly  -  f(x).y*-  fM,(x).ly3  .  fww(x)'yy  =0. 

1!  2!  3!        T.  4! 

or,    f (x)  -  flizhz\  V^iSht¥  ■  0. 
21  4! 

f'(x)  -  fwl(x)y*  ■  0* 
1!  3! 

These  give  us  the  two  equations 

UVx3*1  -7x  +  10)  -  (6x^+3x   -5)y^yy-  o. 

(4x3+3xx-10x  -7)  -  (4x-M)y^=  0, 

Eliminating  y  from  these  two  equations, we  have  after  simplifying 

16x  +24x*-28xV-38x3-32xl~-11x-6  m  0. 

The  roots  of  this  equation  are  -2,  2,  and  four  imaginary  roots. 

2 

Substituting  these  values  of  x  in  the  second  equation  in  x  and  y, 

we  find  for  x  ■  -2,  y  ■  ~£\ .Therefore  the  complex  roots  of  the 

equation  are  -2  +  1,  and  -2  -  1.  For  x  =  £,y  becomes+Ji, which  is 

2  "2 

excluded  since  y  must  be  real. If  we  wished  to  permit  y  to  become 

imaginary, our  roots  would  be  3  *  J  -  1  ■  1  #and  2  -  J.i*i  =  2 

2      2  2      2  2  2 

These  are  the  real  roots  of  the  given  equation* 

k«Method  of  Intersection  of  Curves  in  the  JPlane  of  the  Complex 
Variable  ♦  (Rectangular  Coordinates  , )  *  V/  6 1/ 1  bj  2  °' 

This  method  has  for  it's  main  purpose  the  determination  of 
a  first  approximation  to  the  complex  roots  of  an  equation. The 
knowledge  of  such  a  first  approximation  is  of  great  importance 
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In  the  methods  of  Newton, Schroder, Horner, and  Weddle,as  we  shall 
see  in  Chapter  IV • 

In  the  given  equation 

t(z)  m  0, 

make  a  substitution  as  in  the  method  of  the  preceding  paragraph, 

z  =  x-fiy, 

thus  obtaining  the  two  equations  in  x  any  y 

u  m  0,  v  »  0, 

where  u  represents  those  terms  of  f(x-fiy)  free  of  l,v  those 
involving  i.A  value  of  z  =  x  +  iy  which  is  a  root  of  f(z)  =  0 
must  satisfy  both  equations  u  =  0  and  v  =  0  simultaneously; that 
is, the  graphs  of  these  two  equations  must  Intersect  at  that  point* 

In  order  to  plot  these  graphs, proceed  as  follows: Substitute 
any  constant  value  c,  for  x  in  the  function  u  obtaining  a  function 
in  y  alone •Substitute  for  y  successive  positive  integers  d,,d2,... 
and  determine  the  sign  of  the  function  for  each  value  assigned. A 
change  of  sign  between  two  successive  integers  will  indicate  that 
the  curve  u  =  0  crosses  the  line  x  =  c,  between  them.Repeat  the 
process, using  a  new  value  x  ■  c^.It  is  well  to  use  successive 
Integers  both  positive  and  negative  for  c, ,cit . . . .The  intervals 
on  the  lines  x  =  c,,x  =  cz,...  should  be  plotted  and  approximate 
curves  drawn  through  the  positive  and  negative  ends  of  the 
intervals. Since  y  enters  only  to  even  powers, the  curve  will  be 
symmetric  with  respect  to  the  x-axis.The  true  curve  will  not 
necessarily  lie  between  these  two  curves, but  from  the  nature  of 

in  general  ^quires 

the  plot ,willAlie  wholly  within  the  ME— 1  bounded  by  tlm  llnme 
x  a  cc    and  §  through  which  the  curves  pass. 
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Having  plotted  approximately  the  curve  u  =  0,we  should  treat 
the  curve  v  =  0  in  a  similar  manner. Again  we  have  symmetry  with 
respect  to  the  x-axls,for  y  enters  only  to  the  odd  powers. This  can  be 
readily  seen  if  we  write  u  =  yf(x,yx)  .The  (p-functlon  is  symmetric 
with  respect  to  the  x-axls. Hence  u  =  0  will  be  symmetric  also, since 
the  sign  of  the  factor  y  cannot  affect  the  vanishing  of  <p  .For  all 
values  where  u  ^  0, however, the  function  u  will  have  opposite  signs, 
for  y  taken  positive, from  the  corresponding  values  for  y  negative. 

The  root  of  f(z)  =  0  will  lie  in  the  squares  through  which 
the  curves  of  both  u  =  0  and  v  =  0  may  pass. If  a  smaller  area 
is  desired, it  may  be  obtained  by  using  fractional  values  for  x  and 
y.ln  this  case  only  the  area  included  in  the  squares  in  which  the 
root  lies  needs  to  be  examined. 

Scheffler  fo3has  developed  a  similar  method  for  finding  the 
intersections  of  u  =  0  and  v  =  O.He  substitutes  values  of  x  and 
y  aa  described  above, but  instead  of  determining  merely  the  signs 
of  the  functions  u  and  v,he  computes  their  values  at  each  point. 
These  he  plots, using  the  point  (c,<,d^)as  origin, the  value  of 
u  as  the  abscissa, that  of  v  as  the  ordinate. In  this  way  he  obtains 
curves  traced  by  the  end- points  of  u  and  v,and  thus  finds  more 
accurately  the  points  where  each  curve  u  =  0  and  v  =  0  crosses 
the  line  x  =  oK  . 

This  method  requires  special  care  in  the  case  of  equal  or 
nearly  equal  roots. If  two  roots  were  equal, the  curves  W,  »  0 
and  v  =  0  would  not  cross  each  other  at  the  point, but  would 
merely  become  tangent. In  this  case  one  looking  for  roots  might 
overlook  the  two. One  is  led  to  suspect  the  presence  of  roots 
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by  the  near  approach  of  the  curves  to  each  other; in  this  case 
the  interval  in  which  this  occurred  should  be  investigated  very 
carefully, or  the  equation  examined  for  equal  roots  by  one  of  the 
well-established  methods* In  fact, it  is  usually  preferable  to  get 
rid  of  any  rational  or  multiple  roots  before  attempting  to  find 
complex  roots. 

This  method  is  usually  much  simplified  by  the  knowledge  of 
upper  and  lower  limits  to  the  roots •Knowing  these  limits, one  can 
confine  himself  to  the  intervals  in  which  the  roots  lie  and  often 
avoid  much  tedious  and  fruitless  labor  in  the  parts  of  the  plane 
lying  outside  of  these  boundaries. 

Example:  z*i  zy  -5Z2"  -8z  +10  =  0. 
u  =       x3  -(6y%  5)x^  -(3yx+8)X  +  yv+5yVlO  =  0, 
v  =  (4x3+3x*"  -4xy2-10x  -y*-8)y  =  0. 
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These  Intervals  as  given  in  the  tables  are  plotted  on  the 
following  page, the  area  between  the  two  v-curves  being  shown  in 
black, that  between  the  ti-curves  in  red.Since  the  two  intersect 
in  the  areas  bounded  by 

x  =  -1 ,  x  =  -3,  y  =  1,  y  =  2, 
and  x  =  -1,  x  =  -3,  y  =  -1,  y  =  -3, 

we  know  that  a  conplex  root  of  f(z)  lies  in  each  of  these  intervals. 
Also, since  the  x-axis  is  part  of  the  v-curve,the  figure  shows 
that  real  roots  lie  between  0  and  1 ,and  between  2  and  3. 
The  approximations 

*rr  ly,  =-2+1-1.5 

x-iy,  =-2  -1-1.5 

will  have  real  parts  correct  to  within  1  unit, and  imaginary  parts 
correct  to  within  0.5  unit. 

The  solution  of  this  equation  by  Horner's  method, given  on 
page  41, finds  the  roots  to  be  approximately 
x  +  iy  «  -2.009fl5)+  1. 08 l(99|i, 
x  -iy  =  -2.009(15) -1.08 1(9S)1. 
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CHAPTER  III, 

i 

METHODS  OP  REPEATED  LIMITS.  (  POLAR  COORDINATES  ). 

7. Bernoulli's  Method.  Uj 

Given  an  equation  f(x)  =  0  with  roots  x , ,x2 ,x3, . . .xn  , where 
|x,  |  rjXxl  2|*3j3 ♦  ••^Wf.tlWM  we  have 

m  t       .  ...  T     '#>  6f  (if  -+(if 

(1)    L   ss  l    x,  —   =  x, 


and 


X,  -f  X^  ^  X3  f.M-fX^ 


(i:fef(i;y:-(i:j|i;f  _ 


Build  up  the  equations  whose  roots  are  the  1st, 2nd, 3rd, , , ,  ,l$th 
powers  of  the  roots  of  f(x)  *  0  and  from  the(m-1)th  and(m)th 
equations  calculate  x    by  substituting  in  the  left-hand  side  of  (1) 
the  values  of  the  sums  of  the  (m-1)th  and  (m)th  powers  of  the  roots 
of  the  given  equation .These  sums  can  be  easily  obtained  by  taking 
the  negative  ratio  of  the  coefficient  of  the  second  term  to  the 
coefficient  of  the  first  term  in  each  of  the  respective  equations $ 

A  more  powerful  method     for  determining  these  successive 
sums  of  powers  of  roots  is  given  by  the  use  of  Newton's  Formulae* 
Prom  these  we  have; 

s,  =  -a,,  =  -a, 

bx  m  a,x  -2aif  =  -a,s,-2ai 

s3  =  -a^3a,at-3ai,    ■  -a^-a^s,  -3a  3 
 -•t»A-i.**A*A-3*  ..•-a„_i8,-(n-1)a„_,  . 
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where  bj    «    x^Xj-f        .  •  • .  +  xj|  ,and  where  a , . .  .a„  are  the 


coefficients  of  the  terms  of  f  (x)  =  0,when  a  0  ,the  coefficient  of 
the  first  term  is  1 .For  values  of  s n  ,sA4),..«#  we  have  the  formula 

By  means  of  these  formulae  we  can  calculate  smfor  values  of  m 
as  large  as  we  desire • 

In  case  the  equation  f(x)  =  0  has  two  or  more  real  roots  of 
equal  largest  absolute  value  the  ratio    s  ^    may  not  give  the 
value  of  one  of  the  roots* If  ,say ,of  the  /M^roots  of  equal, 


largest  absolute  value .A of  thuw  are  positive  and  y  negative, we 
have 


«   £££  |x 

SJ-S  -t  1 

•  X*4  1 

s    x 

/**■  ~+  y  1 

and 


The  ratio  a*-*/     would  not  approach  a  limit  therefore  ♦However, 

•——————■»    ■  L 


L  --------    =    L  -------  *  x ^  s  x^~ 88  •••  =  x^v, 


In  case  M-  m  V  #  &   —    becomes  Indeterminate, but  L  — 

still  gives  the  desired  values. 

If  the  two  roots  of  largest  absolute  value  are  conjugate 
eomplex.say/f^6*1^^  and /f^0-c^e);it  is  possible  to  determine 
first  the  angle  0  .  We  may  write  ; 

0m  *  2/^*0  ■+  (x*-*-  fr?-*   *x*  ) 

where       is  the  absolute  value  of  x3, which  is  less  than  f  *,and 
c  ^  is  a  constant  less  than  1 •  The  ratio    s»w/    will  not  approach 


s 


fc 
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a  limit; for  we  have 


Here  the  cosine  functions  in  both  numerator  and  denominator, 
separated  by  the  phase   0  at  all  times, can  each  take  values 
small  enough  to  make  the  tern  In  which  they  appear  too  small 
for  the  second  term  to  be  neglected, no  matter  how  large  m  be 
taken .Moreover  since  each  oscillates  between  1  and  -1,both  num- 
erator and  denominator  will  take  positive  and  negative  values, 
which  may  for  certain  values  of  Q  cause  the  ratio  to  have  different 
signs  and  therefore  not  approach  a  definite  limit. 

If  the  sequence  s* 
( 3 )  sgn  cos  0  , sgn  cos  2  6  , . . . •  sgn  cos  n  $ 

is  known, we  can  approximate  to  6  in  the  following  manner: 

Let  V,  be  the  number  of  positive  terms  before  the  first 
negative  term  of  the  sequence  ;let  Yz  be  the  number  of  successive 
negative  terms  before  the  next  positive  term, and  so  forth. Then, 
since  we  may  assume  without  loss  of  generality  that  6<  Tf  ,we  know. 


from  which 


U1    fl  =  L    UA-or    _  ,  _ — ill — 

More  briefly, if  w^is  the  number  of  changes  of  sign  in  the  first 
m  terms  of  the  sequence  (3), we  have 
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The  sequence 

sgn  s,  ,sgn  s^,  sgn  s„,  

may  agree  with  the  sequence  (3)  in  such  a  manner  that  6  can  be 

f 

determined  from  it  by  the  use  of  (4). If  however  -y-  is  not  suffi- 

J3 

ciently  small, this  is  not  possible. Netto  (  s   p.267)  gives  the 
conditions  for  this  as  being  dependent  on  f% ,  /^and  om  , which  are 
unknown; so  that  the  conditions  are  valueless  in  practical  use, 
Netto  attributes  the  following  modification  to  Jacobi,who 
followed  a  suggestion  by  Fourier: 

liven  an  equation  which  has  A.  roots  of  equal  largest  absolute 
values  with  the  absolute- values  of  the  sums  of  the  (k)th, (k -M )th, . . . 
powers  of  the  roots  x, ,x2 ,x3 , . . .x  sufficiently  large  in  comparison 
with  the  sums  of  corresponding  powers  of  the  remaining  roots  so  that 
we  may  neglect  powers  of  those  roots  smaller  than  x^c^and  write: 

3y^=    Xy    +  X£  -f-  •  #  »  •  ¥  "Xy^ 


We  know  that  the  unknown  coefficients  of  the  equation 
(5)  0 

with  roots  x,  ,xx  ,  ...  x^  can  be  determined  from  the  relations 
used  above, namely: 

•hW*K  ^S/^35  0 


Eliminating  ^  ;Yl"  Y/^  >  from  ^^fc*011  (5)  and  the  set  of  equations 
numbered  (6), we  have 


D9  r  *i 
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x^  '       ...  1 

3/<+/*   8  Hiv-~-l    •     •     •     3  /< 

=  0. 

For 2, as  is  the  case  when  we  have  two  conjugate  complex  roots, 
this  determinant  reduces  to  the  following: 

Laisant  tfu  gives  a  sequence  for  determining  roots  of  equations 
which  is  more  general  than  the  sequence  made  up  of  the  sums  of  the 
successive  powers  of  the  roots. He  uses  Newton's  formula? 

•*m  -h  +  a^S/n  _x -/•  ...  -f  afJs^,_h  =  Op 

and  proves  that  if  a  sequence  has  for  a  law  of  recurrence 

u^4a/uM.ifaiu^i  +  ...  ■+  a^u,.,,  =  0, 
we  may  choose  any  arbitrary  values  ttp,_,  ,U/n_2.f  ..      u^_n  ,and  the 
sequence  built  from  these  can  be  used  for  determining  the  roots  of 
the  equation  with  coefficients  i  ,a, paLf  .  .  •  an,ln  the  same 
manner  as  the  sequence  made  up  of  the  sums  of  successive  powers 
of  the  roots. 

Example:  Find  the  complex  roots  of  the  equation 

xJ  -  3  x*+x-H5  ■  0. 
By  the  use  of  Newton's  formulae  we  get  the  following  sequence: 
3,     3,    7,    3,-13,-77,-  233,  -557,  -1053,  -1431,  -455,  5331, 
23603,    67753,     153001,  273235,   

Selecting  the  last  four  terms  as  aM    ,BAt,  ,s/x,^,and  BMi.3  respectively, 
and  substituting  in  the  formula  (7)  above, we  obtain  the  following 
equation: 
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9792X1-  38729x  +  47690  =  0. 
The  roots  of  this  equation  are    1.977+  0.9861,  1.977  -  0.9861. 
The  correct  values  of  the  roots  of  the  original  equation  are 
2  4-  1,  2  -  1.  This  is  a  very  special  equation, as  the  real  root 
is  1, which  makes  the  sums  of  powers  of  the  complex  roots  grow 
large  very  rapidly  in  comparison.  In  general  this  is  not  the  case, 
so  that  a  larger  number  of  terms  of  the  sequence  will  need  to  be 
found  to  obtain  equally  accurate  results. 

8.3raffe'a  Method.  '< »  '»       3 '>  **>  7 1 ' 

The  method  proposed  by  D.Bernoulli  gave*  only  the  roots  of 
greatest  and  least  absolute  value .In  1837  C.H.Graff e3'developed  a 
method  which  used  many  of  the  same  principles  involved  in 
Bernoulli* s  method, but  which  gave  all  the  roots  of  an  equation 
at  once, with  a  more  rapid  approach  to  the  result. His  method  is  as 
follows: 

Given  the  equation  f(x)  =  0  with  roots  x,#Xit«.««  x„.  The 
product  f(x).f(-i)  ■  0  will  be  an  equation  of  the  (atn)th  degree 
with  roots  x,  ,-x,  ,x2,-Xi,. , . . .  x^-x^.In  performing  the  multipli- 
cation all  odd  powers  of  x  will  disappear, and  the  substitution 

x*  *  x^ 

will  therefore  give  an  equation  of  the  (n)th  degree  in  x' , 

(x«-  xf)(x»  -  if*1  -  x,v)  =  0, 

whose  roots  x/",x^~,  ...  x^    are  the  squares  of  the  roots  of  the 

original  equation* 

The  new  equation  can  be  operated  on  as  was  the  original 

equation, and  the  process  continued  indefinitely, giving  equations 
♦Meyer  extended  it  to  all  roots  later.Math.Annalen,v.33, 1889. p. 51 1 . 
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with  roots 


*  t  2_  #X  ^   ,     ♦  »  .  •    X  ^ 

V*" 


*  i    »X  ^  ,  X  ^  >  X  ^ 

Graff e  then  proves  that  if  all  the  roots  are  real^  unequal  so  that 

Jx,|>  |xi|>|x3J>....>|xJ 
the  equation  of  the  (r)th  powers  of  the  roots  takes  the  form 
( approximately) , 

(1)       x    -  x,x    -+  x,x^-x     -x,XjX3»x     +  •  •••  +  x,x^x3..«x^  =  0, 
for  all  values  of  r  sufficiently  large »His  proof  is  based  on  the 
symmetric  functions  of  the  roots  of  an  equation  whose  coefficient 
of  the  term  of  highest  degree  in  x  is  unity ♦ 

Coefficient  of  x*"'  =  -(x^+x^  +  x^V. .  J  x^  ) 

Coefficient  of  x*^  =  • 

Coefficient  of  x"^3  =  -(x^x^x^xfx^x^  . .  .-^x^x'Jf • .  *4x£a£x«) 

Coefficient  of  x°     ■  {•if  xfxfx3\...  X% 
In  each  case  the  first  term, involving  powers  of  roots  large  in 
comparison  to  those  of  the  following  terms, would  for  Increasing 
r  become  so  large  that  the  remaining  terms  could  be  neglected, 
thus  giving  equation  (1)  above. 

From  equation  (1)  any  real  root  can  be  found .For  example, 
if  x  ^  is  desired, it  can  be  determined  from  the  following  ratio: 
A        coefficient  of  x  xixjX3  •  ••x^fxjJ_ 

J<        coefficient  of  x*1^'  x^x^xj1  ...x£, 

Now  consider  the  case  where  the  two  roots  of  largest  absolute 
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value  are  conjugate  complex. Since  they  have  equal  absolute  value, 
the  assumption  previously  made  that  the  absolute  value  of  x,ls 
greater  than  the  absolute  value  of       no  longer  holds. We  may 
assume  the  following, however: 

I*J>  |x3j  > .....  >  J**}  • 

Equation  ( 1 )  now  becomes  : 

(2)  x    -(x|-fx1)x  -+  XjX^x     -x,XiX3x    -/-  ...»  ±x(x3_...x/1  =  0, 

In  this  case  the  value  of  x,  can  be  determined*  from  the  relation 

jxn,\  -  |x*|  «|V50C| 

In  case  the  complex  roots  are  not  of  greatest  absolute  value, 

the  necessary  assumption  becomes: 

I  *  i  J  >  K  /?•'••>  /«*  | 88  I  *  •••?/»*/• 
Equation  (1)  then  becomes  as  follows: 

i   -±x,x2  ....x^x„  ■  0. 

The  absolute  value  of  x  t<  and  xk+f  can  be  obtained  from  the 
coefficients  of  *****  ,xW<,and  x*~*~'  . 

To  determine  the  argument  of  the  roots  x  (  and  x^  ,we  make 
use  of  the  fact  that  ^    xj  ■  -a,    (assuming  that  a0  =  I.)  If  we 
determine  first  the  absolute  values  of  all  the  roots, we  can  find 
the  real  roots  by  substitution  in  f(x)  =  0  to  determine  the  sign 
of  each. Then, representing  the  two  complex  roots  by  yj(oos  0, -/- 1  sin#, ) 
andf^cosfy-i  sin#,  ),we  have  the  relation: 

(3)  2  /{cos  &,  +x3  -f  Xy7....:f  x„  =  -a, , 

••a  i  •x^^.Xy".  •  •  •  •  _x 
or,  cos  B,  a   -  ^  

 2  n   ,  

*/xi/  and    )x  Jcannot  be  determined  from  the  coefficient  (xj'+x^  ) 
since  Jx,  |-/-|Xi    may  be  greater  than    Jx?-f      | . 
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the.  espon^in^  tc» 

If  the  equation  has  two  pairs  of  complex  roots?equationA{3) 
will  contain  two  unknowns  cos     and  cos  04« In  this  case  the  follow- 
ing additional  condition  is  needed: 

2  i  eo&Q,  j  2  ;  cosft^t   _j_  t-f*  ^  %  a  -an 

This  will  give  both  values  of  $  when  combined  with  (3)  which 
now  becomes  as  follows: 

2^008^^- 2/^0086)^X5--/-  x^  -f-x^  =  -a, 

If  more  than  two  pairs  of  complex  roots  occur, the  following 

fro 

method  will  determine  the  arguments  of  the  roots: 

Let  the  degree  of  the  equation  be  even, equal  to  2m. (If  it  is 

odd, multiply  by  x  throughout, making  it  even)* If  the  equation  is 

divided  by  xM,it  will  take  the  following  form: 
a^x-j^x      ,.i+a2sx  =0. 

If  now  feL6  be  substituted  for  x  and  the  real  and  imaginary 

parts  each  equated  to  zero, the  two  equations  become  as  follows: 
a£cos  mdi  a!cos(m-1  )&-h  . ..  -f  a**,  =0 
b'0sin  m$  +  b*8in(m-1)^  +  . -f  bj,sin#  =  0, 

where       a^  m  a<,/ +a^  m  bj 

•  •••••• 

ai.(*  a^f+U  ft**,  f'1  =K-, 

From  Trigonometry  we  get  the  following  formulae: 
cos  20  =  2  cos*©  -1 , 
cos  30  ■  4  cos30  -3  cos 6 

if  x 

cos  4^  =  8  cos  6  -8  cos#  -f  1 , 

cos  58  =  16  cos5^  -20  cosWs  cos# 


•<1  • 
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sin  28   ■    2  cos 9 
sin  6 

_sin  3  0  =    4  Cosx6)  -  1 . 
sin  a 

sin  4  a  =    8  cosJ0  -  4  cos  0 
sin  0 

sin  5  0  ■  16  cosV0  -12  eosS  4  I. 
sin  9 

These  values  substituted  in  the  two  equations  above  will  give 
for  each  absolute  value  p  of  a  root  two  equations  in  sin  0  and 
cosg     from  which   9  can  be  calculated » (as  in  ptr.HjP.Sj     fi»*\*i  tu 

For  the  case  in  which  more  than  two  roots  have  equal  absolute 
value  a  method  corresponding  to  the  case  of  two  roots  of  equal 
absolute  value  has  been  given. C*Runge  *°  gave  in  1900  the  follow- 
ing method  which  deals  with  a  more  general  case: 

Given  an  equation  f(x)  5  %0X%mrX*  ^  aix""V  ••••^a>,  ■  0. 
If  the  roots  of  this  equation  are  such  that  x,,x^,....  are 
large  in  absolute  value  in  comparison  with        »a$*fi  *•       *h  » 
the  given  equation  may  be  broken  up  into  two  equations: 

a^x^ajX^-'-/-  a^x^V  . .  • -/  a^=  0, 

such  that  the  first  equation  will  have  the  roots  xj9xLt..»  x^, 
approximately, and  the  latter  the  roots  x^,  ,x/^,...«xM  , approx- 
imately. If , in  turn, the  roots  of  these  equations  can  be  divided 
into  sets  one  of  which  contains  only  roots  large  in  eomparison 
with  the  roots  of  the  other  set, the  two  equations  can  again  be 
broken  up, and  the  process  continued  until  the  resulting  equations 
become  of  as  low  a  degree  as  is  desired, or  until  all  the  roots 
of  any  one  equation  become  of  absolute  value  such  that  no  one  is 
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sufficiently  large  In  comparison  with  the  others  to  permit  of 
further  breaking  up  of  the  equation^ i.e. ,all  roots  have  equalA value. 

Conversely, if  each  of  the  coefficients  a6,a,,...       is  large 
in  comparison  with  each  of  *^         >•••  aMtne  equation 

f(x)  g  afl  x V 9l,  xn~ '+  a2xA~V  MM  4  a  n  =  0 
can  be  broken  up  into  two  equations 

a„  x^+ a,  x^V  a^x^"^  a3x/M~^..«  4  a^  •  0, 

a*    x    4-a^,  x  ^   -f  ....4  an  =  0, 
the  first  of  which  will  give  approximately  the  ^  roots  of  largest 
absolute  value  of  f (x)  ■  0, while  the  second  will  give  the  remaining 
h         roots . 

This  gives  a  means  of  simplifying  the  solution  of  any  equation 
with  several  real  or  complex  roots  of  approximately  equal  absolute 
value; for  if  the  equation  of  the  powers  of  the  roots  of  the 
given  equation  is  built  up  as  previously  outlined, the  equation 
will  separate  into  equations  of  lower  degree  which  can  be  solved 
more  readily  than  the  original. 

To  simplify  the  method  of  obtaining  equations  each  with  roots 
equal  to  the  square  of  the  roots  of  the  preceding  equation, we 
have  the  following: 

f  (x)  a  a^x  *  a,  x     y...4a^,=  0 

f(-x)=  a0xA-a, x*~' 4  ...  ±an  =  0. 

Then  f(x)«f(-x)  =  a^xn  -(a^-2a0az)x'li,^V(a^-2a/a37t  2a^av)x  +   » 

+  M)      ( a,>2a^ 2a^aA -2a«, a^j. . . * 2a0 a^)x  ±  

or,     b^x*  4-b/xt    -f  b^x*    -4- bn  =  0, 
where     b0=  a^ 

 b,°  -(alL«2aJaz,)  


4    .     |  • 


\ 


»    |  | 


I 
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(a^-2a/a3+2a0ay) 

b3=-(a^-2  axa/2  a,a^-2  a0a6) 

b  s 

j  H 

(-1)  (a^-2^a/^  2a^^^a^v-2ay^b a/^# •       [2a4  aj 

(-D  < 

With  these 

relations 

the  step  from  one  equation  to  the  next  can 

easily  be  made* 

Example:  find  all  the 

roote  of  the  equation: 

v,  3 
X  -/X 

-5x  -73ff  10  =  0, 

Power 

6f 

roots •  a  o 

• 

&  < 

a  a.                a  a 

1  1 

1 

1  w 

a>1* 

»>   1  " 

a>25          i>  49 

a*  =100 

a  ^=10 

a,a3=14  aiav=100 

a  a,  =20 

2'  1 

11 

59  149 

100 

1 

121 

3481               222^  * 

100^  * 

-118 

jA ( O  1 

200   

2X  1 

3 

403  104^ 

IOO2. 

1 

9 

1623  1C8fc, 

62,  a, 

2~ 

-797 

•  1203  100i 

I00t 

*  222^    signifies  222 

•to1;  I06 2  ,  1o6-1o^  j  1o8^ 

,108-10^  ;etc. 

.  «  *  I 


1  * 


4 
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Power 

of 
roots 

T  h 

a3 

ay 

1 

6353 
-240^ 

144  ^ 
1594  j 

100 /v 

0  w 

100/V 

■ 

t 

174.0  , 

1 00  ,cj 

100  , 

1 

1569 
-348^ 

303  A0 
-79*„ 

°- 

100  jo 

100ic 

1 

i 

— 4  no  . 

100 

1  WW  30 

100  0 

1 

36820 
-448  ^  o 

511  v.. 
38 

100 
0 

100  ^ 

•  * 

I 

1no  / 

100,  , 

1 

64  v^ 

3014  hL 

 0  yO 

100, 

1°0/to 

t 

i 

1 00 , , 

10O  ,  w 

1 

912/7V 

100ai-y 

1 00  x  s  y 

To  carry  the  work  beyond  the  2    power  of  the 

0  ny 
roots  is 

useless, 

as  four 

of  the  coefficients  are  not  affected  by  further  trans- 

f orma t ion* , whll e 

the  second 

coefficient  will 

continue 

to  oscillate 

lndefinltely4Using  the  equation  of  the  2  7  power  of  the 

roots, 

we  have  the  following: 

OS 
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x*   -*1f34Vay  -+3416fCx'  -y- lOO^x1  -+  100,2(,=  0.    x»  ■  x 


*3» 


I/Too7T 


7 


■--  ■  1.9996 

3016 


if* 

By  substitution  in  f (x)  ■  0  we  find  that  both  these  roots  are 
positive . 

/,  «K«  YUMfU    =  2.2356 
x^x^x^x*,  ■  1  -f  1.9996  4  2(2.2358)cos6  =  -1. 

cos  6  =  -  -2i28S§~ 
4.4716 

e   =  153°  26'  tp". 
x,  «    2.2358(cos  153°26,20w4- i  sin  153*26' 20" 


■  -1.9998  +  0.99981 
x^=  -1.9998  -0.99981 


(The  correct  values  of  the  roots 
are  1,  2,  -2  + i,  -2  -  i.) 


9.Hayashi's  Method. 


T.Hayashi  ^  has  given  a  method  for  the  determination  of 
complex  roots, which  he  applies  to  an  equation  of  the  third  degree; 
but  the  method  is  readily  extended, as  he  remarks, to  «ay  equations 
of  any  degree. For  the  cubic, the  method  is  as  follows: 

Given  the  cubic  equation 

a<>x  +  a,xVajX-*  a3«  0, 
with  real  coefficient*; and  let  the  two  complex  roots  be  feci)  an*  fet&. 
If  we  write  t  f  or  eL°  and  b0,b,  ,bx,b3  ,for  &0f>3,&lf2~  .a^y9  9&3 
respectively , we  have  the  following: 
b^Vb,£'"-+bJL£-*b3  a  0, 
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Multiplying  the  latter  by  £Jwe  obtain 

bjlVb^+b^+bo  =  0 
Therefore  £  is  a  common  root  of 
b^x^+bjX^-f  b^x+b3=  0, 
and  b^xVb^x^bjX+boJB  0. 

Combining  to  get  rid  of  the  terms  in  xJ  ,we  have 

(b^-b^Jx**  (bj^-b^jx  ^(bj-b,1  )  «  0, 
of  which   £.    is  a  root.But  C  is  also  a  root .Therefore, 
(t)      (b,b3-bib0)£^  (b^b^bji -4  (b32~-b6v  )  =  0, 

(bjb^-b.bJ^V  d»^b,b0)^  (b3L  -b^  )  =  0. 
Multiplying  the  latter  by  , 
(2)      (*s  •**  )  iV(b^b3-byb0)e-f  (b,b3.bzb0  )  •  0. 
Subtracting  (1)  from  (2), and  factoring, 

B»Z        Mb.bj-b^bjjt^-l)  =0 
Since  for  complex  roots  6  ^  0  and  therefore  (f^  t)  ^  0,we  have 

(b^-b^  )-(b(b3-bJ>0)  *  0. 
Replacing  the  letters  b^b^b^bj  by  their  respective  values, 
this  equation  becomes 

a«f  -a^^4a,a3p  -a3  ■  0. 
This  is  a  cubic  equation  in  p2"  and  can  be  solved  for  real, positive 
roots, among  which  will  be  found  the  absolute  value      P     of  the 
roots    f&  ,  fjL 

The  method  of  procedure  in  the  general  case  from  this  point 
on  is  not  clearly  indicated. No  means  is  given  for  determining 
which  real  positive  root  of  the  last  equation  will  be  the  absolute 
value  of  a  fcoot  of  the  original  equation. In  the  equation  of  the 
third  degree  this  equation/Ais  the  equation  whose  roots  are  equal 
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to  the  products  of  pairs  of  roots  of  the  original  equation, which 
can  have  only  one  real  root , namely, the  product  of  the  two  complex 
roots  of  the  given  equation .For  the  equation  of  the  fourth  degree, 
however,  the  equation  in  ^becomes  of  degree  8;  for  the  fifth  degree 
equation, of  degree  20;for  the  (n)th  degree  equation, of  degree  n.2fl~3; 
so  that  obviously  the  roots  are  no  longer  the  product  of  pairs  of 
roots  of  the  original  equation, No  statement  is  made  as  to  the 
possibility  of  positive  roots  entering  the  equation  in  which 
do  not  give  absolute  values  of  complex  roots  of  the  given  equation* 

Another  disadvantage  of  the  method  is  the  rapid  increase  of 
the  degree  of  the  equation  in  f  for  increasing  n.For  example, to 
solve  an  equation  of  the  sixth  degree  it  is  necessary  to  find  the 
real  positive  roots  of  an  equation  of  degree  48. In  this  respect 
the  method  offers  little  advantage  over  the  methods  of  Lagrange  * 
and  Taylor  ^;for  while  the  equation  to  be  solved  for  real  roots 
can  be  more  readily  built  up  than  by  those  two  methods, it  is  in 
general  of  so  much  higher  degree  that  any  advantage  of  simplicity 
of  formation  is  overcome  by  the  disadvantage  of  a  longer  and  more 
difficult  solution. 
Example;      x^-hx*  -5x^7x  -MO  P  0. 

We  treat  first  the  general  equation  of  the  fourth  degree  as  in  the 

case  of  the  cubic  discussed  above.Placing  ec6/  =  C  and  bd,b,,b^, 

to       3  ]~ 

b3, by., equal  to  j>ra.0,  f  a,  ,  f  az,  f  &j  ,    a^  respectively, and  reducing 
to  an  equation  of  the  second  degree  as  before  we  get  the  equation 
(bxby  -b,b2bj:-  b>J>^b*bz  -b  ,b3by  ^b^b.,  -^b0b>v  -b'b,^  ) 
+  (b3b^  -bob.b^ -b}b3bw4bob,  -b ^b^ +b0b,bxbv  4b0b^b3bv  -b^bj 
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Eliminating  the  second  term  ts§  liuf mm  and  replacing  bj'  by  it's 


value  in  terms  of  a;  and  /  we  obtain  the  following  equation: 

t*Z  )/  -(a03aJf/4-(aJa,a3)//V(aXav  "a^>  h  -*S»J>  f'° 
4(2apafa3av-2aja^)/'  +  (a„a2a5   -a^ajav  -a^)/  ^-{a  ,a3a  v  )/ 

-ta^JfVa^  =  0. 
In  the  given  equations 0  =  1,a,  =  1  ,a  ^=  -5>a  3  =  -7, a  ^  =  10. 
Making  this  substitution  we  get  from  the  equation  above: 

5f'^  -7f'Z-\09f/0- 340/^-1090/^-700/^5000/^10000  =  0. 
The  real  positive  roots  of  this  equation  are  )[5~  and   |/2.There  are 
two  conjugate  complex  roots  of  absolute  value       ♦  The  arguments 
of  these  roots  can  be  found  from  the  quadratic  equation  in  £ 
given  above.  There  are  no  complex  roots  with  absolute  value 


10. Solution  by  Means  of  the  Bquatlon  of  Products  of  Pairs  of  Roots. 


The  following  method  was  suggested  by  the  method  of  Hayashi 
when  applied  to  equations  of  the  third  degree .Hayashi  forms  an 
equation  from  the  given  equation  by  means  of  the  process  described 
in  the  preceding  article. In  the  case  that  the  given  equation  is  a 
cubic, the  new  equation  turns  out  to  be  that  equation  whose  roots 
are  the  products  of  the  roots  of  the  original  equation  taken 
together  in  pairs. For  higher  degrees  this  is  not  the  case;but 
the  lower  degree  equation  suggested  that  there  might  be  some 
value  in  a  general  method  which  gave  such  results  for  all  degrees. 

Thus, given  an  equation  of  degree  n  with  roots  x,....  xn,form 
by  means  of  symmetric  functions  of  the  roots  the  equation  whose 
roots  are 

XjX^f  X(X^,    ....  X,  X  ^  ^ 
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X^X  3  , 


X^X  K, 


This  equation  will  be  of  degree    n(n-1 )    •  Since  the  product  of 

2 

two  conjugate  complex  quantities  is  always  a  real  positive  quantity, 
we  need  only  solve  this  new  equation  for    real  positive  roots. 
Prom  these  can  be  found  the  absolute  values  of  each  of  the  complex 
roots  of  the  original  equation.For  example, if  x^  and  x    are  conju- 
gate complex, we  have  the  relation? 

|x,  |  ■  p\  )fx^xx\ 

For  the  equation  of  the  fourth  degree 
&0xH-+  a;x3-f  a^xS  a3x  -t-  a^  ■  0 
with  roots  x, ,xi,xi,xv,the  new  equation  becomes  as  follows: 

a  0x*^  -a^aax'  ±  a^(a(a3-av)x,^-(a^av-2a<7axa /a^apx* 
+  apa^a ,a3-a¥)xl  -a^a^x'-f-  av    s  0 
whose  roots  are  x,xi,x,x3,x)xv,xix3,xixv,xiXv 

This  method  has  one  disadvantage  of  the  method  of  Hayashi  in 
that  each  positive  root  of  the  new  equation  may  not  give  the 
absolute  value  of  a  pair  of  complex  roots  of  the  original  equation. 
It  has  the  advantage  of  giving  an  equation  of  lower  degree  to 
be  solved  for  real  roots, however jwhile  the  knowledge  that  the  roots 
of  the  new  equation  are  going  to  be  the  products  of  pairs  of  roots 
of  the  original  equation  may  make  it  possible  to  eliminate  at  least 
in  some  cases  the  disadvantage  mentioned  above . 
Example: Find  the  complex  roots  of  the  equation 

x^+Sx^Tx2"  -3x  -10  =  0. 
The  equation  of  the  products  of  pairs  of  roots  is  the  following: 
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x'b-7x'*  ^x^-tldx*3  +50x,;L -700x*  -1000  ■  0. 

The  roots  of  this  equation  are  5»-2,and  four  imaginary  roots .Since 
four  are  imaginary , the  original  equation  must  have  at  least  two 
complex  roots.  Since  5  is  the  only  real  positive  root, it  cannot 
have  more  than  one  pair  of  conjugate  complex  roots. Therefore  the 
original  equation  has  a  pair  of  complex  roots  and  two  real  roots. 
The  absolute  value  of  each  of  the  complex  roots  is  /fT • 

1 1 « Method  of  Intersection  of  Curves  in  the  Plane  of  the  Complex 
Variable.  (Polar  Coordinates)/^  S  *'  7s>  z  0t 

The  method  as  outlined  in  paragraph  6  can  also  be  carried 
out  by  means  of  polar  coordinates .Given  the  equation 

=  adx'V  •/*"~-f  a«  *  0 

make  the  substitution 

x  m  f  (cos  &  -f  1  sin  6  ) 
and  separate  the  resulting  equation  as  before  into  two  equations 
u  g  uff^  )  z  a0f  cos  n&f  a,  f     cos(n-1  )&    •  *««  +%n  =  0, 
t  =  v(f^)  =        sin  nf  -h&,f     sin(n-1  )0-f  sin  e  =0. 

Plot  the  two  equations  u  =  0  and  v  =  0,as  in  the  case  of  rectangular 
coordinates, by  assigning  values  to  f   and  &  and  determining 
intervals  in  which  the  curves  lie. The  roots  of  f(x)  =  0  will  be 
found  on  the  intersection  of  the  two  curves. 

The  chief  advantage  of  this  method  is  that  one  can  determine 
asymptotes  to  the  curves  and  use  these  as  starting  points  for 
plotting. These  curves  have  been  examined  carefully  **'>d'so  that 
methods  of  plotting  can  be  greatly  simplified  by  results  previously 
determined. One  of  these, for  example, is  the  fact  that  the  asymptotes 
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of  either  of  the  curves  u  =  0  or  v  =  0  pass  through  a  common  point 
which  is  the  origin  when  the  original  equation  is  transformed  so 
that  the  coefficient  of  the  (n-t)th  power  of  the  variable  is  zero; 
and  that  they  divide  the  plane  into  n  equal  sectors .We  would 
expect  to  find  a  large  value  of  P  along  a  radius  vector  near  one  of 
the  asymptotes; while  if  we  found  that  there  were  no  real  point 
of  the  curve  on  two  vectors  lying  in  the  same  sector, we  would  know 
that  the  curve  would  not  occur  at  any  point  lying  between  them. 

M.Naraniengar^  in  a  resume  of  methods  of  finding  complex 
roots  of  equations  discusses  this  method  and  gives  plots  of  the 
u  and  v  curves  for  several  equations. 


CHAPTER  IV. 
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METHODS  OP  DOUBLE  LIMITS. 


12. Newton's  Method.  <^  H^,^ 


If  there  is  given  an  approximation  ^  to  a  root  f  of  the 
equation 

t  M  «  o, 

Newton's  method  consists  of  finding  successive  numbers  yx  }}f \y  Y^ ••■ 
each  after  some  one  YH  being  nearer  to  the  root  y  than  the  preceding 
one  of  the  sequence, and  the  limit  of  the  sequence  as  a  becomes 
infinite  being  equal  to  the  root .The  method  is  as  follows: 

If  y la  the  root  and  y  is  the  first  approximation  to  it, 
designate  by  y,  the  difference 

Y'-Y-Y, 

By  means  of  Taylor's  theorem  we  have 

ny^y?)  f  Hh+mw;'+  siLmrf"+  °. 

1 !  2! 

If  yjla  sufficiently  small, we  may  neglect  powers  of  it  above  the 
first ; and  solving  the  equation  above  for  K,  we  find  : 

y/m . 

Then  Yi     Yi    s  Yl  wil1  De  a  ndW  approximation  to  y9  and  can  be 
used  in  determining  a  third  value     ,and  the  process  continued 
until  sufficiently  good  approximations  are  obtained. 
For  complex  rootsvlet  the  approximation  be 
Y,  ■  ftieoaB^l  sin*,)  g  . 
Substituting  in  the  formula  given  above, 


n  -   s  .   •   > 

Vtyi  )  na^'  +  (n-1)a, + a,, 


we  get  the  following  value  of  y,  : 


aa^f/'^coafn-D^-fi  sin(n-1K)-f  + 


Inserting  the  numerical  values  of  the  sines  and  cosines, performing 
the  indicated  multiplications, collecting  terms, and  dividing, we 


One  disadvantage  of  the  method  is  that  it  is  long  and  tedious, 
as  the  convergence  is  often  very  slow ♦ Another  more  serious  one  is 
that  while  the  method  gives  a  sequence  of  values  approaching  the 
root  as  a  limit, it  gives  at  no  time  a  knowledge  of  how  close  to 
the  root  the  approximation  really  is, In  the  case  of  real  roots  this 
difficulty  has  been  overcome  for  a  large  number  of  cases  by 
the  union  of  Newton's  method  with  the  Regula  Falsi  method.By  this 
means  two  sequences  of  points  are  obtained  which  except  for  certain 
special  cases  approach  the  root  from  opposite  sides, at  least  after 
a  certain  number  of  approximations  have  been  made; and  thus  is 
obtained  an  interval  between  the  points  of  the  two  sequences  in 
which  the  root  must  lie. No  such  improvement  seems  to  have  been 
made  for  the  case  of  complex  roots. To  do  so  would  probably  be 
rather  difficult, if  possible  at  all, since  the  roots  may  lie  any- 
where in  a  plane  instead  of  on  a  known  line* 

The  conditions  for  a  given  approximation  leading  to  any 
particular  root  have  been  carefully  worked  out  for  the  case  of 


find 


to 
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of  real  roots  and  Cay ley  ,v  has  done  the  same  for  the  complex  roots 
of  the  quadratic  equation .Beyond  this  little  seems  to  have  been 
done  In  regard  to  finding  regions  in  which  approximations  must 
lie  to  lead  to  a  particular  complex  rootAor  in  fact  to  lead  to 
any  one  of  the  complex  roots  of  an  equation. 
Example:    Find  the  complex  root  of  the  equation 

f(x)  =  x3  +  x^-x  -M5  =  0, 

o  o 

with  a  given  approximation  «  3(cos  60-hl  sin  60  )  =1.5+"2.6i. 

f '(x)  =  3x+2x  -  U 

+15 

/      27(cos  180^-1  sin  l8cf)+9(cos  120+  1  sin  120°)  -3(cos  60°iBln66) 

v;-  -  -  - 

27(cos  120+1  sin  120° ) -+6(cos  60  +  1  sin  60° )  -  1 
-27  -  4. 5  +  7. 81  -  1.5  -2.6i  -f  15 
-13.5+-23.41  +  3  +  5.21  -  1 
=  -0.374  -  0.4791 

Yi     ■  vi+  y,'  «  1.5  +2,6i  -  0.374  -  0.479 

=  1.126  -+-  2.1211  ■  2.402(cos  62V  -+  1  sin  62°  2' ) 

Use  for  the  new  approximation  ^=  2.40(cos  62tf2'  -f  i  sin  62°2/). 

/  (2.4oftcos  186*6+1  sin  186V)  +  (2.4oftcos  124*4  +  1  sin  124%'  ) 
JTx*  -   -(2.40) (cos  62*2^1  sin  62g2  )-f  15 

3(2.4oftcos  124*  4  +  1  3ln  124° 4)+  2(2.40)(cos  62'2+i  sin  62°2)-1 
=  -0.113  -  0.109 
^3=  1.126 -t  2.121  1  .  0.113  -  0.109  i  z.  1.013+-2.012  i. 
The  true  roots  are  1+  21, and  1-21. 

13. Horner's  Method.  Ly 
Given  an  equation 

fix)  m  0, 

with  a  root  y  -  ^  +  i(3  ♦Assume  that  the  first  figures  o^and  ^<3(  of 
♦.See, however, Schroder, ^  P.353. 
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the  real  and  imaginary  parts  respectively  of  this  root  are  known* 
Diminish  the  roots  of  the  equation  by  the  amount  ^  -+  l  p  .This  can 
be  done  by  synthetic  division  as  in  the  case  of  real  roots*  The 
coefficients  will  become  complex  however. The  first  figures  of  the 
root  of  the  new  equation  will  be  equal  to  the  second  figures  of  the 
root  of  the  original  equation. Let  them  be^and  /3i#The  roots  of 
the  new  equation  should  now  be  diminished  by  o£ ^  ftand  the  process 
continued  until  a  sufficient  degree  of  approximation  is  reached. 

When  the  root  of  one  of  the  transformed  equations  becomes  of 
absolute  value  less  than  1 ,the  first  figure  of  the  root  can  be  fauuwty 
obtained  very  readily  from  the  last  two  terms  of  the  equation. 
Assume  this  equation  to  be 

(a^ib'Jx'V  -f(aj_^  ib*_,  )x'-t(a\-hlb\)  =  0. 

Since  the  absolute  value  of  the  root  is  less  than  1 ,we  neglect 
powers  of  x1  above  the  first, and  solve  for  x'  from  the  last  two 
terms , obtaining 

If  this  approximation  does  not  give  the  exact  first  figure  of 
the  root, the  error  will  be  corrected  in  the  process  of  finding 
the  next  succeeding  figure  of  the  root* 

Geometrically  the  method  is  equivalent  to  the  following: 
Assume  a  root  of  the  equation  at  the  point  P  in  the  plane  of  the 
complex  variable. Let  the  point  P.  represent  the  point  whose 

P0i.*\t     wt\0SC    coord.  L.  H\Al^^  -he- 

coordinates  are  the  first  figures  of  the  root  P;P^the^ first  two 
figures ; etc .The  successive  transformations  are  equivalent  to  the 
transformation  of  the  origin  to  the  points  P,,Pit...  successively, 


* 


0 
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the  axes  being  translated  parallel  to  themselves ♦The  roots  of  the 
successive  equations  become  smaller  and  smaller  in  both  their  real 
and  their  imaginary  parts ;and  the  root  of  the  original  equation 
is  given  by  the  sum  of  all  the  real  parts, plus  the  sum  of  all  the 
imaginary  parts, of  the  roots  of  the  transformed  equations* 


Example:  Find  the  complex  roots  of  the  equation 

x^/x3  -5xr  -8x  +  10  =  0* 
From  paragraph  6  we  find  the  first  figures  of  th#  roots  to  be 
-2  4i,and  -2-i.  We  proceed  to  find  the  former  root  by  diminishing 
the  roots  of  the  equation  by  -24- i.  —  y{  . 

1  f-1  .5  -8  -MO  \  -2+i 

-8 


-2    -hi  -tt 


-31 


4  11 


4-  2i 


-i 


1-1     -+1  -4 
-2    4i  4-4 


-3i 
-71 


43 
4-10 


4  21 
f-20i 


1  -3 
-2 


2i  e 

41  +7 


-10i 

-111 


1  -5  m 

-2  41 


-t  2 


-i 


i  13 


4221 
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-211 


1  -7  +41 

The  first  transformed  equation, therefore, is  as  follows: 
x,-(7-4l)x,V(7-2ll)x,X4'(13  4-22i)x*4-2-i  ■  0. 
We  find  x1  from  the  last  two  terms ♦ 


at  m.z  • 


41. 

-2  H*  i 

Yx  -x*  m   s  0.00  +  0.081 

13  +  221 

We  diminish  the  roots  of  the  equation  by  this  amount* 

1  -7-Mi    +7      ^211      413        +  221       +2        -1         lo.QO  0.081 

0  +  0.081+0.3264-0.561  +  1 .7248  +0.53391- 1  ♦ 8027+1 .17801 

1  -7    4.081  6.6736-21 .561+14.7348+22. 5339i|+0. 1973+0. 17801 

0  +0.081-0.3328-0.561  +1 .7696  +0.50721 

1-7  +4. 161+6. 3408-22. V2i +16.4944+23.041 U 

Without  carrying  the  work  farther, we  can  obtain  a  third  approximation 

from  the  last  two  terms. 

-0.1973-0.17801  m     v  .  . 

y  -  X"  *   m  -0.009(15)  +0.001(99)1 

16.4944+23.04111 

The  root  of  the  original  equation, therefore, is  the  sum  of  the 

approximations  ^  +  JV+X  ■  -2  +  1  +0.00-^0.081  -0.009(15^0.001(99)1 

=  -2.009(15)  ~h 1.081(99)1 

H.Weddle's  Method. 

The  method  of  Newton, as  described  in  article  12, assumes  the 
knowledge  of  an  approximation  to  a  root  of  an  equation, and  seeks 
to  obtain  a  better  approximation  by  addition  to  the  known  one; that 
is, if  yt  is  the  known  approximation, the  value        <f/(lf  Y/^s 
properly  defined)  will  be  a  better  one. For  an  infinite  number  of 
operations, the  method  develops  into  the  infinite  series 

.  U I + r,'*aL  t  v:  •..  •• + in  t ;.  •  •  -  r. 

Thomas  Weddle  has  developed  a  somewhat  similar  method .Given  an 
approximation  ]f  to  a  root  \  ,he  seeks  a  value      such  that  the 
product  X(,_f  ^")wlll  be  a  better  approximation  to   ft  than  is  Y,  ♦ 
For  an  infinite  number  of  operations, he  gets  the  infinite  product 
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He  proceeds  as  follows:  Let  ^  be  a  root  of  the  equation 
f(x)  ■  0,and  assume  that  Jfj  is  an  approximation  sufficiently  near 
to  j^.If  we  transform  the  equation  by  means  of  the  substitution 
(D  x  =  ^x' 

we  will  obtain  a  new  equation  P(x)  ■  0  whose  root  ^  *  will 

01 

not  be  much  different  in  absolute  value  from  unity .This  root  may  be 

written  Jf=  (1  +  Y"), where  "   y'Srill  be  a  very  small  quantity"  in 

absolute  value* Expanding  by  means  of  Taylor's  theorem, we  get 

f(i-v  n  ■  F(1)+r,(1))'V  r,(\)y"\  .  .  .  =0. 
//  2! 

Neglecting  powers  of  %  above  the  first .we  find  an  approximation  to  I  , 

Y"  F(lJ 
f  m  

F'(1) 

An  approximation  to  the  root  of  the  transformed  equation  , therefore, 

(/    nCi*J  appKi**- 'fn  Co  the 

is  (1+^  ).TheAroot  of  the  original  equation, from  ( 1 )  ,1s    Xk  (  1  ~*~  //  ) 

The  method  is  satisfactory  in  certain  cases ♦Diffleul ties 
enter, how ever, in  regard  to  the  region  in  which  an  approximation 
must  lie  in  order  to  lead  to  a  root  of  the  equation .Moreover, if 
\Y'I  is  approximately  1 ,then  Y  " 9  ( defined  by  the  relation  Ym  1  +  ?  ') » 
is  not  necessarily  very  small  in  absolute  value. In  fact, for  the 
case  where    y/  =  i,  (and  therefore   |  y'\  =1),   J  )"  isl/i"  .Other 
examples  are  easily  found.These  difficulties  need  to  be  taken  care 
of  before  the  method  can  be  safely  used. 

"Example: Given  f(x)  =  x3+xx-x  -f  15  ■  0,and  x  «  1 .5  +  2.6i  an 
approximation  to  ft  complex  root .Find  the  root. 

iy=  1.5  +  2.6i  =  3(cos  60°-M  sin  6o" )  approximately. 
Substituting  x  =  3(cos  60 V  i  sin  60*)XjWe  get  the  equation 

27(cos  18oV  i  sin  180°)x'-h  9(cos  120*Vi  sin  120° )x' 

 -3(cos  60° -hi  3 in  60°)x'  -f- 15  -  0 

♦This  example  is  given  to  explain  the  method  of  procedure  rather 
than  to  show  the  faults  of  the  method. The  ^glutlpn  is  satisfactory 
in  this  special  case. In  general, the  resultsSA«lar£> uncertain. 
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Replacing  the  cosines  and  sines  by  their  lmmorloal  values  and 

dividing  both  sides  of  the  equation  by  -3, we  get  the  new  equation 

F,(x0=9x,3+(1.5  -  2.6i)x,:L-t(0.5  +  0.87i)x!  -5  =  0, 

.//    •„  Ml)         91-1.5-2.61  +  0.5+0.871-5 

iftjU-ft--  -  -  ------  —  

Fj(1)  27+3  -5.21-/- 0.5 +0.871 

■  iO. 815 +0.0331 

(1  +o(,+  c(3()  a  0.815  +  0.0331    g  0.815(cos  2°19  +  1  sin  2°l9  ) 
Substituting  x'=  0.815(cos  2°19  -fi  sin  2°  19  )xH  we  get  the  equation 
27(0. 815)^(008  186*57  +1  sin  186° 57'  )xwJ~ 

+9(0. 815)*~(cos  124*38+1  sin  124°38')xf,i 

-3(0.815)(cos  62dl9/"fi  sin  62°19/)xVl5  ■  0. 
Replacing  the  cosines  and  sines  by  their  numorloal  values  and 
simplifying, we  get  the  equation 

(4.835  +1  0.589)xw3-(1. 136-1. 638i)xw  +  (0. 379+ 0.724i)xH  -5  $  0. 
tMirf-  ^ftx)  ■  -  — —  *  -0.0776  +  0.02251. 

(l-Hfj+'A)  *  0.9224+0.02251  «  0.9224(cos  1°24Vi  sin  \° 24  ) . 
The  root  of  the  original  equation  is  given  by  the  product 

■  3(cos  60 Vi  sin  60°) (0.815) (cos  2°  19+1  sin  2°  19') 

.(0.9224)  (cos  1°24+i  sin  1°24' ) . 
=2.255(cos  63°41+i  sin  63°4l') 
=  0.9997  +  2.0211 
The  true  root  is  1  i-  21. 

15. Method  of  Infinite  Series.  '     '  ' 

E«McClintock  H<f  has  made  use  of  the  Laplace  series  v3in 
developing  a  "method  for  finding  all  the  roots  of  an  algebraic 
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equation  simultaneously**  .For  the  trinomial  equation 

x   *  w  +nax 

he  derives  from  the  Laplace  series  the  following  series: 

where  y  is  a  root  of  the  trinomial  equation, and  u/  is  any  one  of 
the  (njbh  roots  of  w  • 

He  gives  the  following  condition  as  being  necessary  and 
sufficient  for  convergence: 

|a  |<  k    (n-k)  |      for  n  always  positive. 

If  for  any  given  equation  the  series  does  not  converge, he  considers 
the  equation  in  two  forms .First, dividing  by  na  and  rearranging 
terms, he  gets  the  equation 

x     ss  -  — -  -j-  

na  na 

or,placing  —  w    =  vut  and    1    ■  (n-k)a,, 
na  na 

xM^(  -Mn-k)a;x1 

Second, dividing  the  given  equation  by  x"~*,he  gets  the  equation 
x    ■  na  +  lo  x 

He  shows  that  if  the  original  equation  does  not  give  a  convergent 
series, each  of  these  two  equations  will  do  so, the  first  one 
giving  n-k  roots, the  second  the  remaining  k  roots, of  the  original 
equation ♦ 

For  the  n-term  equation  he  uses  the  more  general  series 


where  the  given  equation  was 


x   cui  -fna^>(x)       (in  general,  a  =  1). 
Having  given  an  equation  for  solution, his  first  step  is  to 


For  other  methods  of  solving  trinomials, see  Weber, p. 397. 


* 


» 


t.ij  -.7 
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examine  it's  fitness; that  is, the  convergency  of  the  series  derived 
from  it .A  criterion  for  testing  this"fltness"is  not  given, it  being 
suggested  that  it  can  be  done  by  inspection. If  the  equation  is  not 
"fit", a  linear  transformation  exists  (so  the  author  states  without 
proof)  that  will  render  it  "fit" .This  transformation  should  be  found 
and  applied. 

Let  us  assume  that  the  equation  after  having  been  made  "fit"  is 
the  following: 

a^x-fa;x    +  p •  m-f-m^Z  -H    •••-fa/x  a^x  -fa„=  0. 

The  next  step  is  to  select  "dominant"coefficlents.This  is  done  in 
accordance  with  the  conditions 

I.  a^,a^,are  always"dominants" 

I I.  "Dominants "must  be  large  in  comparison  with  other  coefficients. 

III.  Three  coefficients  ay<,a^  ,aM, (coefficients  of  terms  in  x  " ~f<  9x"  s 
x^'respectively)  are  called  Success ive"dominants "if  |aT~"/ia 
decidedly"  larger  than  |a^_/<  a^"jt  That  is, this  relation  must  hold 
for  any  set  of  three  succes8lve"dominantsf., 

The  "dominants"  of  the  equation  having  been  found  to  be  a09a^9 
a^9a^,9  ...  a^,the  equation  is  replaced  by  a  set  of  equations 
obtained  as  follows: 

Equation  1:  Divide  the  given  equation  by  atfx"  'and  rearrange 
the  terms  into  the  form  x  H  =  u/+ka^>,(x)  • 

Equation  2:  Divide  the  given  equation  by  ax<x  ~       and  rearrange 
the  terms  into  the  form  xi_/<  =      K  ■  (i  -k)a^(x). 

Equation  3: Divide  the  original  equation  by  a^x  'and 
rearrange  the  terms  into  the  form  x"     =  u/~(m-i)a<2>(x). 


In  general, the  f -functions  will  have  negative  exponents , and  a  may  be 
taken  equal  to  1 .  ;  


'  i' 


1 
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The  series/yabove  Is  then  applied  to  each  equation.All  the  roots 
can  be  obtained  from  these  equations .The  first  equation  will  give 
a  different  root  for  each  of  the  (k)th  roots  of  urK  , giving  k  roots 
in  all;the  second  equation  will  give  1  -  k  roots ;and  so  forth, the 
entire  set  giving  therefore  all  the  roots  of  the  original  equation* 

The  method  seems  to  be  of  considerable  practical  value* 
licClintock  solves  several  problems , taking  successive  ones  from  a 
textbook1* on  the  Theory  of  Equations  to  avoid  criticism  of  having 
selected  special  types  of  examples. If  the  theoretical  foundations 
could  be  strengthened  by  means  of  criteria  for  fitness  of  equations, 
some  certain  way  of  selecting  suitable  transformations , or  a  more 
exact  method  of  determining  domenant  coefficients, the  work  would 
be  of  much  more  value  and  usefulness. 

P.Lambert  HljliX  has  developed  a  similar  method  which  also  makes 
use  of  the  Laplace  Series. For  the  three  term  equation  he  proceeds 
as  follows  I 

Let  the  given  equation  be 


Form  the  three  equations 

( 1)  a^xV  a^x'y-l-a, 

(2)  aflxVaw<x%^ 


0, 


(3)  a0x*Vfa,  an 


0. 


These  may  be  written  as  follows: 


(3«)  %  uL  ys — -x  / 
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Expand  each  of  these  functions  on  the  right  into  a  power  series  in  y 
by  means  of  the  Laplace  series, which  is  as  follows? 

If*''*3     x  =  fO+y<p(x)] 
Then  P(x)  =  F(f(z))-hy(p(f(z))^|F(f(z)))  -|j^[fUTzT)^(f U)$]t. . 

4  [?lFU))n^(f(z)e]  + ... 

where  P(x)  represents  a  root  of  the  equation  f(x)=0  and  F(f(z))  is 
a  root  of  the  equation 

adxV  a^  ■  0, 
the  roots  of  which  are  assumed  to  be  known. 

If  after  expansion  y  is  made  unity, the  resulting  P(x)  becomes 
one  of  the  roots  of  the  given  equation, provided  the  series  is 
convergent .To  each  of  the  three  forms  above  corresponds  a  series. 
For  each  of  the  n  roots  of  the  binomial  equation  a0xA-/-aM  =  0,one 
of  these  series  will  be  convergent. Each  root  will  give  a  different 
root  of  the  original  equation, so  that  of  all  n  roots  are  used, all 
the  roots  of  the  orignial  equation  will  be  obtained. 

The  method  works  out  very  well  for  the  three  term  equation, 
and  the  conditions  for  the  convergence  of  the  series  are  simply 
stated. For  equations  with  more  than  three  terms , however , the  process 
becomes  very  much  involved. The  solution, for  example, of  a  six- term 
equation  depends  upon  the  solution  of  a  five-term, a  four-term, a 
three  term, and  a  two-term  equation.The  solution  of  the  four-term 
equation  demands  the  convergence  of  a  doubly  infinite  series, in 

aj1  Sent*. 

which  every  term  of  an  infinite  series  isAinf initeaA  five-term 
equation  demands  the  convergence  of  a  triply  infinite  series, and 
so  on. 

While  theoretically  the  method  may  be  of  some  value, the 
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difficulties  presented  would  seem  to  be  of  such  a  nature  as  to 
render  it  almost  useless  in  practical  work  for  determining  the 
roots  of  equations  of  more  than  three  or  four  terms* 

16. Schroder's  Method •  UH 

E.Schroder  UH  made  an  attempt  to  generalize  the  methods  of 
Itui'atlus  in  order  to  bring   many       methods  into  one  group, He 
proceeds  as  follows: 

Given  an  equation 

f  (x)  =  0 

and  an  approximation  °< ,  to  a  root  o(  •  The  provlem  is  to  find  a 
function  F  such  that 

will  be  a  better  approximation  to  o{  than  4,  .  Then 

43=  F(<)  =  P(F(ot,  )) 
will  be  better  than  either  <& ,  or  4x9 and  in  general 

H^it-i  )  =  F(F(....(F(4()  )) 

will  be  better  than  any  preceding   <4H  .The  further  condition  is 
required  that 

L         m  <  • 

It  we  have  a  function  that  satisfies  these  conditions, so 
that, at  least  from  3ome  point  on, each  successive  application  of 
F  gives  us  a  closer  approximation  to  the  root  < ,we  have  an 
algorithm  for  the  solution  of  equations  for  both  real  and  complex 
roots. 

One  example  of  an  F  function  is  given  by  Newton's  method: 

d.i  .4  .  t<*,)    =  »(«,) 
f  u ,  ) 
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• 

Each  auccesaive  application  of  P  gives^a  value  more  nearly  equal 

to  the  root  than  any  of  the  preceding  approximation, at  leaat 

from 

aome  point  on. 

The  conditiona  aa  staged  above  which  the  function  F  muat 

aatiafy  can  be  more  compactly  atated  aa  follows : 

I.           >  U  -  «U  >U~<i\ >  ■  •  '  '  >M~*U^  •  •  '  ' 

II.    l_     U-  <AJ  =  0. 

It  is  obvious , also, that  we  muat  have  from  I 

X'.FW)  ■<<  . 

If  we  expand 

F(oUO  =  T«)-f  £F!W+  ItVH)^  •••• 

and  if  we  place  F(<0  =  ck   and  F(  <A  +  £  )  ■       ,we  may  diaregard 

powers  of  t  above  the  first, obtaining  the  relation 

or,  t'F'W*  , 

which  is  equivalent  to  saying  that  if  the  new  correction  4ro( 

is 

to  be  in  absolute  value  less  than  the  original  correction  £ ,the 

following  relation  must  hold: 

II*.  (F*  (ok)  |  <  1. 

Thus, in  the  case  of  Newton1 s  method, where 

FUJ  =  |tiJL 
ru,) 

we  have  (since  f  (  ^  )  =0  ) 

IS    F(  <*)  =  c{  , 
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II'.  F'U)  =  I^Jii  =  I  -  ^y^-f  '      "  °- 
Schroder  also  takes  up  other  algorithms  for  the  solution  of 
equations  which  include  a  number  of  other  methods  as  special  cases. 
One  of  these  is  the  following: 

<*,f("0 

where  A  is  a  parameter  which  can  take  any  value. If  A  becomes  zero, 
we  see  that  the  algorithm  reduces  to  Newton's  method .Another 
algorithm  is 

Schroder  obtains  another  class  of  methods  by  combining  an 
infinite  number  of  operations  into  one  step. These  methods  are  not 
discussed  in  this  thesis  on  account  of  the  difficulty  of  clearly 
presenting  the  matter  in  brief »It  also  seems  that  their  practical 
value  may  be  less  than  their  theoretical  value, although  several 
important  known  methods  (such  as  Newton's  method)  are  obtained  . 
from  the  general  method  by  specialization. 

17 .Additional  Methods. 

A  number  of  other  methods  for  solving  equations  with  complex 
roots  have  been  developed. One  in  particular  that  should  be  mentioned 
is  the  Fiirstenau  method  of  infinite  determlnants.Gunther,33'' 33  , about 
1870, examined  this  method  and  developed  it  to  some  extent .Concerning 
his  first  article,the°Jahrbuch  iiber  die  Fortschritte  der  Mathematik 
saysTMDer  daraus  gezogene  Schluss,dass  der  wahre  Werth  der  Wurzel 
zwischen  zwei  auf elnanderfolgenden  Naherungswerthen  liege, 1st  aber 
nicht  begriindet,da  u.A.der  entsprechende  Satz  nicht  von  allgemeinen 
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Kettenbruchen  gilt."  Later, In  1877,Gunther3*solves  the  equation 

(x-mT  =  x^4xV6xV4x  -M  =  0 

and  finds  a  root  to  be  -  £  approximately , which  he  says  is  a  good 

3 

approximation  since  it  reduces  the  left-hand  side  to  0.012  . 
Naegelbach     has  also  discussed  the  method  and  has  solved  several 
problems, but  the  process  seems  too  complicated  for  practical 
purposes* 

Gonggrijp  *°  has  recently  published  a  method  related  to  Graff  e's 
method  which  seems  to  offer  difficulties  in  the  determination  of 
successive  derivatives. These  appear  to  make  it  impractical  for 
general  use, though  it  may  have  considerable  theoretical  value. 

Naraniengar     gives  a  geometrical  interpretation  to  several 
methods  that  have  been  discussed  in  this  thesis* 

Weber  12~  gives  a  method  in  connection  with  his  proof  of  the 
Fundamental  Theorem  of  Algebra  which  is  of  theoretical  importance 
rather  than  of  practical  value* It  consists  in  finding  circles  in 
which  roots  lie, and  finding  smaller  circles  within  these  which 
contain  roots, and  so  forth. It  is  given  to  show  that  the  Fundamental 
Theorem  proof  leads  to  a  method  of  finding  roots  as  well  as 
proving  their  existence, 

Runge/1  Cohn,/fe  and  others  have  given  methods  based  on  the 

/InalyCCc  of  4  Compel  ^i-'it/e 

Theory  of AFunct ions A which  could  not  be  treated  here. The  problem  of 
separation  of  roots  (that  is, the  determination  of  intervals  or 
areas  in  which  only  one  root  lies)  is  too  extensive  to  be  discussed 
in  this  article.Graphlcal  methods,  ^> b^      , several  electro«?ohemioal 
methods  *1>Wt*,s'3  based  upon  work  by  Lucas, and  a  mechanical  method 
developed  by  L.Torrea  *  r  lie  beyond  the  limits  of  this  thesis* 
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♦Murphy's  Algebraical  Equations  contains  an  explanation  of  the 
work  of  Simpson. 


In  general, only  the  first  page  of  the  discussion  of  interest  to 
this  thesis  is  given. 
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